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ABSTRACT 

The relative alignments of mid-infrared traced Galactic bubbles are compared to the orientation of 
the mean Galactic magnetic field in the disk. The orientations of bubbles in the northern Galactic 
plane were measured and are consistent with random orientations - no preferential alignment with 
respect to the Galactic disk was found. A subsample of Hii region driven Galactic bubbles was 
identified, and as a single population they show random orientations. When this subsample was fur- 
ther divided into subthermal and suprathermal H ii regions, based on hydrogren radio recombination 
linewidths, the subthermal Hii regions showed a marginal deviation from random orientations, but the 
suprathermal H ii regions showed significant alignment with the Galactic plane. The mean orientation 
of the Galactic disk magnetic field was characterized using new near-infrared starlight polarimetry 
and the suprathermal Hii regions were found to preferentially align with the disk magnetic field. If 
suprathermal linewidths are associated with younger H ii regions, then the evolution of young H ii 
regions is significantly affected by the Galactic magnetic field. As Hii regions age, they cease to be 
strongly linked to the Galactic magnetic field, as surrounding density variations come to dominate 
their morphological evolution. From the new observations, the ratios of magnetic-to-ram pressures in 
the expanding ionization fronts were estimated for younger H ii regions. 

Subject headings: Hii regions — Infrared: ISM — ISM: bubbles — ISM: magnetic fields — radio lines: 
ISM — Techniques: polarimetric 



1. INTRODUCTION 

Mid-infrared (MIR) objects called "Galact ic Bub- 
bles" were cataloged bv lChurchwell et ahl (2006') and the 
sky projections of these three-dimensional bubbles were 
shown to be preferentially elliptical. What forces cause 
these bubbles to show non-circular shapes? External 
magnetic fields are one possible explanation, since or- 
dered magnetic fields can apply anisotropic pressure. To 
test this hypothesis, the orientations of a subset of Hii 
region driven Galactic bubbles were compared to predic- 
tions of the mean Galactic magnetic field orientation. 

Anisotropic bubbles can have a number of causes: stel- 
lar motions. Galactic shear, anisotropic driving forces 
(e.g., bipolar outflows), expansion into a non- uniform 
medium, or the presence of a magnetic field, which 
adds anisotropic pressure (see lBis novatvi- Kogan fc Silich 
119951 for a thorough review). Detailed models for the 
expansion of stellar wind-driven bubbles into a uni- 
form interste l lar me dium (ISM) were fi r st dev eloped by 
ICastor et al.l (fT975h and IWeaver et all (flQTl . The ef- 
fects of magnetic fields on bubble expansion, in particu- 
lar how the thickness changes at different points on the 
bubble's exterior shell, were c onsidered by ^Ferrier e et al. 
(i991) for supernovae and iBisnovatvi-K ogan fc Silich 
( 1995 ) for wind-driven bubbles. iTomisakal (|1992i ) quanti- 
fied the effect of a uniform magnetic field on a supernova 
remnant, predicting that the remnant should elongate 
along the magnetic field direction. This prediction has 



been confirmed for supernovae fe.g.. lGaensle3ll998D . but 
n ot for wind-dr i ven or radiation-driven bubbles. 

IWeaver et al.l (|1977D showed that neither stellar mo- 
tions nor Galactic shear can account for the elonga- 
tion of Galactic bubbles. The lifetimes of early-type 
stars responsible for creating these bubbles are much 
shorter than the timescales needed for stellar motion or 
Galactic shear to significantly warp the bubbles. Large 
stellar motions (e.g., runaway O stars) are expected to 
create cone-sha ped bubbles, which were not seen by 
iChurchwell et al. (2006). Planetary nebulae can exhibit 
bipolar nature s, but these a.re likely cau sed by close bi- 
nary svstem s (iBalickl 119871: iSokeii I2006D or stellar mag- 



netic fields (! Jordan et al. I2005f). but not by external 
magnetic fields (jSokeii 120051: iSabin et al.ll2007f ). A high- 
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mass stellar bipolar outflow typically only exists inside 
the s tar's dense natal material (but see lMcGroarty et al.l 
I2OOI for parsec-long flows fro m lower-mass star s) and 
can be readily identifled (e.g., [Chambers et al.l 12009) . 
Based on Herschel Spac e Telescope column density maps, 
[Anderson et al.l (|2012b[ ) found that Galactic bulDble in- 
teriors have low densities. If stellar driven, they have 
evolved beyond the earliest stages of star formation in 
which bipolar outflows typically exist. 

Magnetic fields may affect the evolution of Galactic 
bubbles via the fields providing anisotropic pressure, 
since it is easier for charged particles to move along 
field lines (i.e., lower pressure) than perpendicular to 
field lines (i.e., higher pressure). Only charged species 
sense magnetic fields, but they are coUisionally coupled, 
to varying degrees, to the neutral species. Therefore, the 
presence of ordered magnetic fields (ordered on scales 
larger than the bubble diameter) can allow bulk flows 
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along the field while inhibiting such flows along orthog- 
onal directions. In the presence of an ordered, external 
magnetic field, uniform expansion of a G alactic bubble 
should be warped into an ellipsoid (jTomi saka 199^. 

The outward flow of material making up the bub- 
ble must be driven. Possible mechanisms include 
combin ed stellar winds (IChurchwell et al.ll2006[ ). super- 
novae (iHelfand et al.l l2006D. or expanding H II regions 
(|Anderson et al.ll2011[ ). Since the evolution of bubbles 
likely depends on their energy sources, this work aims to 
avoid such dependence by a narrow focus on H ii region 
driven bubbles. 

Starlight polarimetry, which traces magnetic fields via 
dust alignment, toward each Galactic bubble was drawn 
from early access to th e Galactic Plane Infr ared Polar- 
ization Survey fGPIPS: [Ciemens et all2012cf ). Initial at- 
tempts to directly probe the magnetic fields of individual 
bubbles with GPIPS near-infrared (NIR) starlight po- 
larimetry failed because the bubbles are predominantly 
beyond the distances readily probed by GPIPS. How- 
ever, a statistical analysis of the relative alignment of a 
subsample of Galactic bubbles with the mean Galactic 
magnetic field orientation (which itself is predominantly 
parallel to the Galactic plane) does show evidence for 
correlation, leading to a possible Hii region evolution- 
ary sequence. More turbulent (possibly younger) H ii 
region driven, elongated bubbles are better aligned with 
the Galactic plane than less turbulent (possibly older) 
H II region bubbles. This observationally-driven result 
supports a scenario for the evolution of expanding H ii 
region driven bubbles into an external magnetic field. 

In Section 2, the properties of Hii regions coincident 
with Galactic bubbles are summarized. The method used 
to measure Galactic bubble orientations is presented in 
Section 3 along with the orientations found. In Section 
4, starlight polarimetry from GPIPS is used to estimate 
the mean magnetic field orientation of the Galaxy. Sec- 
tion 5 compares Hii region properties with the relative 
alignment of the elongated bubbles and external mag- 
netic fields. A notional explanation for these observa- 
tions is presented in Section 6 and conclusions are pre- 
sented in Section 7. 

2. H II REGIONS ASSOCIATED WITH GALACTIC 
BUBBLES 

The Gree n Bank Telescope Hii Region Discovery Sur- 
vey fHRDS: lBania" et al. 2010; Anderson et al. 2011) is a 
targeted, 3 cm wavelength hydrogen radio recombination 
line (RRL) and continuum survey of the inner Galac- 
tic plane (343° < i < 67°, \b\ < 1°) that overlaps with 
GPIPS and the Galactic Legacy Infrared Mid-Plane Sur - 
vey Extraordinaire (GLIMPSE; Bcniam inerall 120031 ). 
HRDS targets were selected by searching for overlapping 
24 ^m emission from the Spitzer Spa ce Telescope MIPS 
Galactic Plane Survey (MIPSGAL; ICarev et al.l 120091) 
and coincident 20 cm continuum emission from ei ther the 
NRA O VLA Galactic Plane Survey (VG PS: Stil et al. 
2006^ or NRAO VLA Sky Survey (NVSS; ICondon et a" 
1998). These two target selection wavelengths (24 fj,m 
and 20 cm) were chosen because dust grains in Hii re- 
gions absorb stellar UV photons, which reemit at thermal 
infrared wavelengths, while the plasma surrounding the 
central ionizing star produces free-free thermal emission 
at centimeter wavelengths. HRDS simultaneously mea- 



sured and averaged seven hydrogen recombination line 
profiles (H 86 a to H 93 a, but excluding H 86 a because 
of confusion with higher order RRL transiti ons) towards 
each target (|Balsedl2006l : iBania et al]l2010t ). As of this 
writing, 603 discrete recomination line components from 
448 lines of sight ha ve been detected, an d their radio 
p roperties cataloged (lAnderson et all 1201 It) . 

lAnderson et~aLl (12011) compared the 603 HRD S RRL 
components with the 134 lChurchwell et al.l ()2006[ ) Galac- 
tic bubbles cataloged in the Northern Galactic plane. 
Thirty-three HRDS RRL components with coincident 
with twenty-seven Galactic bubbles (four bubbles have 
two coincident components and one bubble has three co- 
incident components). Thus, the Galactic bubbles are 
assumed to trace the perimeters of H ii regions expanding 
into the interstellar medium. These twenty-seven Galac- 
tic bubbles constitute the subsample which will be stud- 
ied in detail in later sections. 

Selected data from HRDS for the thirty-three RRL 
components coincident with Galactic bubbles are listed 
in Table [T] The observed radial velocity, RRL linewidth, 
and peak RRL fiux towards each of these Galactic bub- 
bles are listed in columns 2, 3, and 4, respectively, in 
the Table. Several sightlines contain multiple RRL com- 
ponents, which are uniquely identified (e.g., N23a and 
N23b) in column 1 of the Table. As seen in the third 
column, several objects show large linewidths. An elec- 
tron temperature of 10'* K results in a hydr ogen thermal 
hnewidth of 22 km s^* (B rown et al.lll978[ ). hence large 
linewidths are suprathermal. 

Using the published HRDS line-of-sight velocity for 
each of the thirty-three recombination line components 
coincident with Galacti c bubbles , kine matic distances 
were calculated using the lClemensI (|T985h rotation curve. 
All Hii regions were assumed to be located at the far 
kinematic distance or the tangent point, where appropri- 
ate, as reasoned below. An intrinsic velocity dispersion 
associated with rando m cloud motions in the Galaxy was 
also included. Though fClemensI (|1985[ ) reports a cloud- 
cloud velocity dispersion of 3 km s~*, th e more conser- 
vative value of 5 km s~* ()Burtonl 119761) was adopted. 
The published uncertainty in the line radial velocity was 
added in quadrature to the random cloud dispersion to 
estimate upper and lower bounds for each kinematic dis- 
tance. This propagated uncertainty was added and sub- 
tracted from each line velocity (RVlsr i (Trvlsr from 
column 2 in Table [1} and used to calculate upper and 
lower uncertainty bounds on each kinematic distance. 
The mean difference between these distance bounds and 
the central kinematic distance was adopted as the kine- 
matic distance uncertainty (column 5 in Table [1]). 

The goal of HRDS was to improve the census of Galac- 
tic Hii regions, especially at large distances. Since the 
nearby population of Hii regions is well known, HRDS 
ignored them. This study was intentionally restricted to 
the HRDS data set because of the uniformity of their 
measurements of the hydrogen RRLs. So, not all Hii 
regions are at the far kinematic distance, but there is a 
very high probability that all the H ii regions taken from 
HRDS arc. 

For sevent een of these t h irty-th ree HRDS RRL 
components, lAnderson et al.l (|2012al ) used Hi emis- 
sion/absorption to resolve the near/far distance ambigu- 
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TABLE 1 

Galactic Bubble H II Region Properties 



Bubble 


RVisr ^ 


Lincwidth^ 


Peak H Line Flux^ 


Dist. 


Luminosity 


a 


b 


Spectral 


Name 


(km s"-"-) 


(km 


s 




(mjy) 


(kpc) 


(Jy kpc2) 


(pc) 


(pc) 


Type 


(1) 


(2) 


(3) 




(4) 




(5) 


(6) 


(7) 


(8) 


(9) 


Nil 


54.0 ± 0.1 


16.5 


± 


0.1 


360 


± 


60 


11.63 ± 1.16 


49 ± 13 


4.57 


3.28 


03 


N20 


39.1 ± 0.2 


19.5 


± 


0.5 


130 


± 


20 


12.72 ± 1.27 


21 ± 5.3 


4.20 


3.16 


06 


N23a 


42.6 ± 0.6 


23.2 


± 


1.3 


390 


± 


30 


12.51 ± 1.25 


61 ± 13 


1.42 


1.31 


03 


N23b 


61.9 ± 7.7 


35.6 


± 


9.9 


390 


± 


30 


11.49 ± 1.15 


51 ± 11 


1.30 


1.20 


03 


N25 


37.8 ± 0.1 


19.9 


± 


0.2 


170 


± 


20 


12.81 ± 1.28 


28 ± 6.5 


3.10 


2.02 


05 


N27a 


60.4 ± 0.1 


18.2 


± 


0.3 


190 


± 


40 


11.54 ± 1.15 


25 ± 7.4 


3.56 


2.95 


06 


N27b 


118 ± 0.3 


22.7 


± 


0.8 


190 


± 


40 


8.14 ± 0.81 


13 ± 3.7 


2.51 


2.08 


06 


N31a 


114.3 ± 0.3 


28.4 


± 


0.8 


340 


± 


80 


7.77 ± 0.78 


21 ± 6.3 


1.61 


1.24 


05 


N31b 


41.9 ± 0.3 


23.6 


± 


0.8 


340 


± 


80 


12.40 ± 1.24 


52 ± 16 


2.56 


1.98 


03 


In 




20.0 


± 


0.2 


170 


± 


10 


Q Qfi -1- n 
c5.yu m u.oy 








07 
\J 1 


N50 


67.7 ± 0.1 


21.7 


± 


0.3 


450 


± 


40 


10.46 ± 1.05 


49 ± 11 


4.84 


4.29 


04 


N53a 


43.3 ± 0.1 


19.6 


± 


0.1 


480 


± 


50 


11.63 ± 1.16 


65 ± 15 


2.50 


2.47 


03 


N53b 


101.6 ± 0.8 


22.6 


± 


1.9 


480 


± 


50 


7.27 ± 0.73 


25 ± 5.8 


1.57 


1.54 


05 


N56 


77.4 ± 0.1 


21.0 


± 


0.3 


270 


± 


20 


9.29 ± 0.93 


23 ± 5.0 


2.72 


2.67 


06 


N57 


30.1 ± 0.5 


22.7 


± 


1.1 


30 lb 


10 


12.24 ± 1.22 


4.4 ± 1.7 


1.10 


0.89 


08 


N60 


50 ± 0.2 


17.3 


± 


0.4 


70 ± 10 


10.84 ± 1.08 


8.2 ± 2.0 


2.02 


1.51 


07 


N62 


62.9 ± 0.1 


18.5 


± 


0.3 


310 ± 40 


9.92 ± 0.99 


31 ± 7.3 


4.09 


3.80 


05 


N66a 


38.1 ± 1.0 


31.0 


± 


2.0 


90 ± 20 


11.36 ± 1.14 


12 ± 3.5 


1.52 


1.22 


06 


N66b 


68.6 ± 1.7 


24.3 


± 


4.2 


90 ± 20 


9.35 ± 0.94 


7.9 ± 2.4 


1.25 


1.01 


07 


N66c 


92.0 ± 0.7 


19.1 


± 


2.1 


90 ± 20 


6.94 ± 0.69 


4.3 ± 1.3 


0.93 


0.75 


BO 


N67 


57.5 ± 0.1 


21.2 


± 


0.1 


470 ± 30 


10.06 ± 1.01 


48 ± 10 


2.20 


2.02 


04 


N73 


60.9 ± 0.2 


28.2 


± 


0.4 


210 ± 20 


9.16 ± 0.92 


18 ± 4.0 


2.06 


1.92 


05 


N75 


42.1 ± 0.2 


21.4 


± 


0.5 


40 ± 10 


10.46 ± 1.05 


4.4 ± 1.4 


1.89 


1.55 


08 


N80 


20.7 ± 0.2 


30.7 


± 


0.5 


620 ± 


80 


11.18 ± 1.12 


77 ± 18 


6.22 


4.52 


03 


N90 


70.5 ± 0.1 


20.8 


± 


0.2 


360 ± 60 


6.14 ± 0.61 


14 ± 3.5 


2.73 


2.71 


07 


N92 


62.5 ± 0.2 


17.0 


± 


0.6 


140 ± 40 


7.29 ± 0.73 


7.4 ± 2.6 


3.59 


2.44 


07 


N95 


52.5 ± 0.1 


20.6 


± 


0.2 


610 ± 30 


8.12 ± 0.81 


40 ± 8.3 


4.97 


3.96 


04 


N96 


-44.4 ± 0.1 


29.3 




0.3 


450 ± 40 


14.65 ± 1.46 


97 ± 21 


1.88 


1.55 


03 


N98 


56.9 ± 0.1 


21.0 


± 


0.2 


390 ± 40 


6.96 ± 0.70 


19 ± 4.2 


2.31 


2.29 


06 


N105 


-1.1 ± 0.2 


18.3 


± 


0.4 


200 ± 30 


10.62 ± 1.06 


23 ± 5.6 


2.83 


1.86 


06 


NllO 


7.9 ± 0.4 


25.7 


± 


0.9 


70 ± 30 


9.59 ± 0.96 


6.4 ± 3.0 


1.60 


1.48 


08 


N115 


23.9 ± 0.3 


29.9 


± 


0.7 


1250 ± 800 


8.14 ± 0.81 


83 ± 56 


7.58 


5.97 


03 


N122 


45.5 ± 0.3 


19.4 


± 


0.6 


110 ± 30 


4.74 ± 0.47 


2.5 ± 0.8 


0.73 


0.51 
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^ Data taken from I Anderson et al.l ||2011I '|. 

ity. All seventeen of the RRL components were either at 
the far kinematic distance or at the tanget point. The H i 
spectra towards the other HRDS RRL components con- 
tained only weak Hi features and no conclusions could 
be drawn. This finding lends support to the assump- 
tion above that all of the thirty-three H ii region driven 
bubbles are at the far kinematic distance. 

With the HRDS RRLs and the kinematic distances to 
the thirty-three Galactic bubbles, additional properties 
can be derived for these H ii region driven bubbles, in- 
cluding the intrinsic luminosities and physical sizes of the 
major and minor axes of the bubbles (columns 6, 7, and 
8 in Tabled. 

The spectral type of the massive star creating each H ii 
region driven bubbles was estimated from the total RRL 
flux and the kinematic distance estimate for each RRL, 
using the Lyman cont i nuum photon em i ssion rates from 
iSteriiberg et al.l (|2003D . lAnderson et al.l ()2011| ) predicted 
the expected HRDS RRL flux density from Hii regions 
hosting stars of these spectral types for a range of dis- 
tances (0 - 23 kpc; shown as their Figure 12), assuming 
that a single star is responsible for all of the ionizing 
flux. From the previously determined peak line inten- 
sity (in mJy) and full- width at half- maximum (FWHM) , 
each RRL component was integrated, assuming a Gaus- 
sian line profile. These spectral type estimates are listed 
in column 9 in Table [TJ 



3. GALACTIC BUBBLE ORIENTATIONS 

iChurchwell et all (2006) fit ellipses to their visually 
identified Galactic bubble sample, as seen in 8 /im 
GLIMPSE images. However, the position angles of the 
bubble major axes ('B-GFA,' in Galactic coordinates) 
were not retained. This key bubble parameter was 
needed for comparison of bubble orientations to the di- 
rection of the Galactic plane, and ultimately to the av- 
erage projected Galactic magnetic field orientation. 

Therefore, an algorithm was developed for fitting B- 
GFAs from GLIMFSE 8 fim images. Its action is il- 
lustrated for the bubble N21 in Figure [TJ In the left 
panel, the 8 /xm emission from GLIMPSE (greyscale) is 
shown, in Galactic coordinates. The goal of the algo- 
rit hm was to fit the inner ellipse parameters reported 
bv IChurchwell et al.l (|2006( ) to the 8 /im emission on the 
inner boundaries of each bubble and so recover the B- 
GPAs. 

First, bright point sources in each image were identi- 
fied with the DAOPHOT FIND routine (,Stetsonil987[) 
and masked, using 18 x 18 arcsec squares. A lSobell (|T97i 
edge enhancement operator was applied to the images. 
This operator calculates the horizontal (G^:) and verti- 
cal {Gy) intensity gradients across each pixel and returns 
an approximation of the amplitude of the total intensity 
gradient (G). In practice, these gradient images are gen- 
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Fig. 1. — (Left panel) Greyscale GLIMPSE 8 /im image of bubble N21 with the final ellipse fit (ellipse in both panels) and the direction 
of the bubble position angle (diagonal dashed line in both panels). (Right panel) Greyscale GLIMPSE 8 fim image of N21 after application 
of a Sobel edge enhancement operator (Sobcl 1978). Points indicate the first 2a datum along each of the 250 radial rays from the bubble's 
center that were used for fitting the ellipse. 

erated by convolving ((g)) the image (A) with two differ- 
ent kernels: 



-1 +1 
-2 +2 
-1 +1 



A; G, 



-1 -2 -1 

+1 +2 +1 



100 - 



(1) 

These two resulting images are combined to create the 
full gradient image G = ^ G^ + Gy. Original image (A) 

regions of constant flux return zero gradients. 

Circular regions centered at the reported bubble cen- 
ters, with radii equal to the reported inner minor axes, 
were masked in each resulting Sobel image (G). The re- 
sulting Sobel image for bubble N21 is shown in the right 
panel of Figure [TJ 

Vectors of pixel values were extracted from the masked 
Sobel image along 250 equally spaced (in azimuthal an- 
gle) radial rays from the reported bubble center. The 
means and standard deviations of values contained in 
each data vector were found. Starting from the center of 
the image and moving radially outward, the location of 
the first datum exhibiting more than a two sigma devi- 
ation above the average value for that vector was iden- 
tified. These 250 locations are shown by the diamond 
symbols in the right panel of Figure [T] These {£, b) loca- 
tions were the basis for new ellipse fitting. 

T he routine MPFIT ELLIPSE from the MPFIT pack- 
age (|Markwardtil2009l ) was used to fit one ellipse to the 
locations of the 250 points described above for each bub- 
ble. The center coordinates, major axis, and minor axi s 
were all fixed and taken from IChurchwell et al.l (j2006[ ). 
leaving B-GPA as the only fit parameter. The result of 
the fit for bubble N21 is shown as the ellipses in both 
panels of Figure [T] The B-GPA is shown as dashed lines 
in the Figure. 

This procedure was applied to all 134 bubbles in the 




0.3 0.5 0.7 0.9 

Bubble Eccentricity 

Fig. 2. — The dependence of fitted Galactic bubble position angle 
uncertainty on bubble eccentricity. The solid line is a power law fit 
to the data. Seven bubbles were fit by hand and their estimated 
uncertainties are shown by the diamonds. Two of the seven bubbles 
had very similar eccentricities and their symbols overlap. 



North ern Galactic plane reported by IChurchwell et al.l 
(|2006[) . The resulting B-GPAs are listed in Tabled along 
with their uncertainties. The uncertainties in B-GPAs 
were found to depend on bubble eccentricity, as revealed 
in Figure [21 For small eccentricities, bubbles shapes and 
position angles are more difficult to ascertain. For bub- 
bles with very small eccentricities, the fitting returned 
B-GPA uncertainities of 180°. The behavior was well-fit 
by a power law {apA = 2.76 ecc"^"^), based on an F- 
test, and the resulting fit is shown by the solid line in 
Figure [2j This same behavior was recovered when only 
the subsample of thirty-three H ii region driven bubbles 
was considered. 

The presence of very bright stars in the GLIMPSE 
field-of-view as well as foreground infrared dark clouds 
caused the fitting procedure to fail for seven bubbles, 
which are identified in Table [3l For these bubbles, B- 
GPAs were fit by hand and the uncertainties (diamonds) 
were conservatively estimated using the published bub- 
ble eccentricities and the power law fit (solid line), with 
an additional 10° of uncertainty. 

The bubble identifiers, eccentricities, and the major 
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axes of the inner bubble wall (all from iChur chwell et alj 
120061) for all 134 northern Galactic bubbles are listed in 
columns 1, 2, and 3 of Table H The B-GPAs and their 
uncertainties are listed in column 4 of that Table. 

4. STARLIGHT POLARIMETRY TOWARDS GALACTIC 
BUBBLES 

All available early access near-infrared starlight po- 
larimetry towards towards Galactic bubbles was ob- 
tained from GPIPS. GPIPS is an H-band (1.6 fim) lin- 
ear imaging polarimetry survey of the northern Galac- 
tic plane (1 8° < £ < 56°, |6| < 1°) using the Mimir 
instrument (jClemens et al.|[2007[ ) on the 1.8 m Perkins 
Telescope outside of Flagstaff, AZ. GPIPS observations 
and data reduc tion pi pelines are discussed in detail in 
iClemens et"all (HoHBIfl). The data used here were ob- 
served between 2007 May and 2011 June. While most 
Galactic bubbles were observed, significant gaps in cov- 
erage surrounding many of them remained at the time of 
this analysis. 

Originally, a goal was to measure stars through the 
interiors of bubbles, through regions exterior to bub- 
bles, and (if possible) through the thick edges of bubbles. 
With this goal in mind, only bubbles with inner major 
axis extents greater than one arcminute were initially 
considered, because of the expected GPIPS polarimet- 
ric sky sampling rate of about one star per square ar- 
cminute. However, this angular sampling rate would not 
provide statist ically sign i ficant samples of stars for most 
of the .Churchwell et all ()2006[ ) Galactic bubbles. Fur- 
thermore, the twenty-seven Galactic bubbles with dis- 
tance estimates are typically beyo nd the GPIPS distanc e 
horizon of approximately 7 kpc ([Clemens "eFall [20123) . 
limiting the utility of starlight polarimetry to probe the 
magnetic fields around those objects. 

Instead, to test the role of magnetic fields in bubble 
asymmetry, bubble orientations were compared with av- 
erage magnetic field directions in the G alactic disk. This 
change in scope allowed study of all 134 lChurchwell et al.l 
12006) bubbles in the Northern Galactic plane instead of 
only those bubbles with the largest angular sizes. 

Lacking direct measures of the magnetic field orien- 
tations surrounding each bubble, the bubbles were in- 
stead assumed to be embedded in a Galactic-scale, or- 
dered magnetic field. The disk of the Milky Way has 
a large-scale, toroidally-dominated magnetic field (e.g., 
[Heiles 2000; Van Eck ct al. 201 1: Pavel et al. 201 2). Al l 
of the Galactic bubbles in the IChurchwell et all ([2006f ) 
sample are within one degree of the Galactic midplane, 
placing them well within the toroidally dominated zone. 
For any disk-symmetric Galactic magnetic field, the pro- 
jection of the toroidal component will exhibit NIR back- 
ground starlight polarization position angles ('P-GPA'), 
in Galactic coordinates, of abou t 90°, for nearly all 
Galactic longitudes and distances ([Pavelll20TT[ ). 

This approximation breaks down for lines-of-sight 
along the local magnetic field direction. In the Solar 
neighborhood, this occurs at Galactic longitudes £ — 
90° -I- p an d £ — 270° + jo, where p is the magnetic 
pitch angle. iPavel et al.l ([20121 ) found p = -6 ± 2°, sim- 
ilar to valu es found by other studies using a variety of 
techniques (IValled [1981 [Han fc Qiaol [l99l IHeilesI [19961: 
[Beca[2007[:[Pshirkov et al.ll2011l) . For this magnetic pitch 
angle, the longitudes of vanishing polarization become 
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Fig. 3. — Histogram of the weighted average GPIPS polariza- 
tion position angles (P-GPA) towards each of the 97 bubbles with 
GPIPS polarimetry. The vertical dashed line at 85.3° shows the 
unweighted average of these measurements. The sharpness of the 
distribution suggests that P-GPA= 85° is representative of average 
orientation of the projected Galactic magnetic field in the disk. 

£ = 84° and 264°, well outside the longitude range of 
GPIPS and where these Galactic bubbles are located. 

W eaker vertical and ra dial magnetic field components 
(e.g.. lFerriere fc Schmitt ''2b00) will perturb the dominant 
toroidal magnetic field. These perturbations will cause 
P-GPA to differ from perfect alig nment with the Galac - 
tic plane (see example deviation in lClemens et al.ll2012aD . 
Therefore, GPIPS starlight polarimetry was used to char- 
acterize the mean P-GPA for the Galactic magnetic field 
within several kpc of the Sun, and this characterization 
was assumed to apply to all locations within the Galaxy. 

Towards each Galactic bubble, all available GPIPS 
starlight polarimetry data within three times the outer 
major ax is radius (Rnut) of each bubble center, as re- 
ported bv lChurchwell et all ([2006f ) , were collected. These 
starlight polarization position angles were combined into 
one weighted mean P-GPA around each bubble, and this 
value is listed in column 8 of Table [3] To keep a few 
very low uncertainty measurements from dominating the 
weighted mean, all starlight polarization position angle 
uncertainties were fioored to 5° or greater. 

The unweighted average starlight polarization <P- 
GPA> towards all of the bubbles, listed in column 8 
in Table [31 is 85.3° with a dispersion of 10°. An F-test 
showed that assuming a constant value for the magnetic 
field orientation across this longitude range was appropri- 
ate and no higher-order longitude-dependent terms were 
needed. Thirty-six Galactic bubbles, N1-N19 and N123- 
N134, are outside of the GPIPS region. Therefore, only 
97 bubbles had GPIPS measurements, and the histogram 
of these <P-GPA> is shown in Fig. [5| 

This mean value agrees with thoretical predictions 
(e.g.. 'Wielebinski & Kraus^ Il993t iFerriere fc SchmitO 
20001: TMoss et al. 2 010) and previous observation s (e.g., 
Mathewson fc Fordl[l970t iHeilcs 2001 iMen efall [20081. 
Analysis of NIR starlight pol arimetry in the outer 
([Pavel et al.l [20T2[ ) and inner ([Clemens et al.l l2012aD 
Galaxy, also supports the conclusion that a single mean 
orientation is a reasonable characterization of the Galac- 
tic magnetic field in the disk. 

5. ANALYSIS 

Using the measured bubble properties from Section [3] 
and the assumption of a constant Galactic magnetic field 
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Fig. 4. — Histogram of the absolute difEerence between the bubble 
major axis orientation (B-GPA) and the average projected Galactic 
magnetic field orientation (P-GPA=85°). 

orientation, the absolute difference between the bubble 
orientation and the average projected Galactic magnetic 
field orientation (AGPA = |B-GPA-P-GPA|) can be 
used to test whether magnetic fields can explain the ob- 
served bubble eccentricities. Two limiting cases are: (1) 
magnetic fields completely determine the bubble orienta- 
tion, or, (2) magnetic fields have no effect on the bubble 
orientation. If the former, then the AGPA histogram 
should be sharply peaked around AGPA = 0°. If mag- 
netic fields have no effect on bubble orientation, the his- 
togram of AGPA should be fiat, consistent with random 
bubble orientations. 

Figure |3] shows the histogram of AGPA for the 134 
Northern Galactic bubbles: the distribution appears to 
be fiat. The cumulative distribution function (CDF) of 
the AGPA data was compared to the CDF of a flat dis- 
tribution, using a KS test, to find a 91% chance that the 
AGPA data were drawn from a fiat distribution. This 
fails to reject the null hypothesis of no preferential align- 
ment between the average Galactic magnetic field orien- 
tations and bubble elongations. When the bubbles are 
considered as a single populations, bubble major axes 
show no preferential alignment with the average Galactic 
magnetic field orientation, and bubble eccentricity is not 
caused by external magnetic fields. AGPA was also ex- 
amined against estimated spectral type of the star pow- 
ering each Hii region, and no correlation (i?^ = 0.0015) 
was found. Yet when subsamples of bubbles were cre- 
ated, correlations emerged, as next described. 

5.1. An H II Region Magnetic Evolutionary Sequence 

Several of the hydrogen RRL components listed in Ta- 
ble [1] are suprathermal suggesting that additional line 
broadening mechansims (e.g., pressure broadening, large 
scale systematic motions, and turbulence) are important. 
At the frequencies of the RRL observations (9 GHz), 
and assuming an electron density of 10^ cm~^ at a tem- 
perature of lO'* K in the line emitting regions, pressure 
broadening shoul d be 1.3% of the thermal broadening 
(|Keto et a l.' 2008, and references therein) , therefore pres- 
sure broadening can be ignored. The effects of large 
scale systematic motions and turbulent motions in the 
gas cannot be disentangled with these observations, so 
their combined effects on the RRL linewidths will be con- 
sidered in the following analysis and collectively refered 
to as turbulence. Assuming a temperature of 10** K in 




Fig. 5. — Measured inner major axis for bubbles along 33 lines- 
of-sight with HRDS recombination line observations as a function 
of hydrogen radio recombination linewidth. The vertical dashed 
line at 22 km s~^ represents the thermal linewidth expected for a 
10'' K gas. The major axes of bubbles N115 and N80, assumed to 
be at their far kinematic distances, were also calculated for their 
near kinematic distances, the two estimates are connected by the 
dotted lines. 



the line emitting regions, the lAnderson et al.l (|2011[ ) RRL 
linewidths can be decomposed into thermal and turbu- 
lent components. 

In Kolmogorov-like turbulence ()Kolmogorovl |1941|) . 
turbulent energy is injected into a system at some char- 
acteristic length scale. The energy is redistributed to 
smaller size scales, reaching a "Kolmogorov" energy 
spectrum of E{k) oc k~^/^, though no assumption is 
made about the actual spectral indicies in the Hii re- 
gions probed here. In the system, energy is dissi- 
pated at some small, resistive, length scale. There is 
a monotonic flow of energy from larger to smaller length 
scales while maintaining the t urbulent energy spectrum. 
[Robertson fc GoldreichI (|2012[ ) considered the effects of 
an expanding or contracting volume containing a tur- 
bulent gas. Their simulations showed that an expanding 
turbulent gas will experience 'adiabatic cooling' and that 
the amplitude of the turbulent energy spectrum will de- 
crease as the volume increases. Under these conditions, 
the turbulent velocity components should decrease with 
time, while the lO'* K thermal linewidth component is 
unchanged. Perhaps suprathermal linewidths are asso- 
ciated with younger, turbulent H ii regions and smaller 
linewidths are associated with older, less turbulent Hii 
regions. 

Figure [5] plots the physical size of the thirty-three H ii 
region driven bubbles (derived from the kinematic dis- 
tance estimates) against the HRDS linewidths of those 
H II regions. A few interesting features are seen in Fig- 
ure [5l The data seem to break into two groups at a 
linewidth of approximately 22 km s^^ (shown by the 
vertical dashed line). At smaller linewidths, a wider 
range of physical Hii region major axes is seen; at larger 
linewidths, the distribution of major axes is generally 
narrower and at smaller physical sizes. Figure [5] shows 
that radius is not an indicator of linewidth. Only phys- 
ically smaller H ii regions (except where there is a dis- 
tance ambiguity; e.g., N80 and N115) exhibit suprather- 
mal linewidths. The full range of physical sizes is repre- 
sented by the H ii regions showing subthermal linewidths. 

Physical size, by itself, should not be a good Hii region 
age indicator because the growth rate of an Hii region 
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will be strongly influenced by the energy output (radi- 
ation and stellar wind) of the cent r al sta r. Model stel- 
lar atmospheres by [Sternberg et al.l ()2003D show that the 
ionizing photon rate can vary by two orders of magnitude 
in going from 0-type to early B-type stars. This causes 
any radius-age relation to become degenerate with the 
spectral type of the central star, in the sense that young 
Hii regions with massive stars will have the same radii 
as older Hii regions driven by less massive stars. While 
an estimate of the spectral type of the star in each of 
H II region was made, this estimate assumed a single ion- 
izing star was driving the Hii region. For the actual 
Hii regions, multiple massive stars could contribute to 
the ionizing fluxes and so would be expected to further 
confuse or mask any age-radius relation. 

To quantify any effect of spectral type on the other 
observed Hii region properties, a number was assigned to 
each spectral type (3 for 03, 4 for 04,. ..,10 for BO). The 
observed hydrogen linewidths and spectral types showed 
a weak linear correlation (i?^ = 0.15), with earlier- type 
stars showing slightly larger linewidths. The correlation 
is too weak for spectral type to account for the break 
seen at a linewidth of 22 km s~^. 

In Figure [3 two outliers (bubbles N80 and N115) 
are seen in the high linewidth region, exhibiting inner 
major axes greater than 6 pc. Both of these bubbles 
were assumed t o be a t their far kinematic distances. 
[Anderson et al.l (|2012al ) used H i spectra to confirm that 
N80 is at the far kinematic distance, but the Hi spectral 
quality was too low to break the distance ambiguity for 
N115. Their physical sizes were recalculated assuming 
they were instead at the near kinematic distance (1.47 
kpc for N80 and 1.96 kpc for N115). If these bubbles are 
at their near kinematic distances, then they agree with 
the overall trend of two groups separated by a character- 
istic linewidth. Nevertheless, the far kinematic distances 
were adopted for the following analysis. 

Since the temperature of an H ii region i s expected to 
be cl ose to 10"* K over most of its lifetime ()Brown et al.l 
119781) , the larger linewidths are likely caused by turbulent 
motions (large scale, systematic motions and small scale 
turbulence) in the line emitting regions of each bubble. 
Also, the physical size of an Hii regions is a monoton- 
ically increasing function with time, which may reach 
an equilibrium size but should never decrease. Together, 
the small physical sizes and suprathermal linewidths sug- 
gest that the H ii regions with larger linewidths may be 
younger. 

Figure [S] shows AGFA as a function of linewidth. 
The vertical dotted line drawn at 22 km s~^ shows the 
linewidth beyond which turbulence must be significant. 
Beyond 22 km s~^, there appears to be better alignment 
(low AGFA) between each bubble and the mean Galactic 
magnetic field. 

Two subsamples were created and analyzed for relative 
alignment of the bubbles and mean magnetic field ori- 
entation; bubbles with subthermal linewidths (Av < 22 
km s~^) and suprathermal linewidths (Av > 22 km s^^). 
For a uniform distribution of relative alignments (AGFA 
independent of linewidth), the average value of AGFA 
should be 45°. The mean weighted AGFA values for 
the subthermal and suprathermal subsamples linewidths 
are 60° ± 8° and 16° ± 6°, respectively. The uncertainty 
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Fig. 6. — Dependence of AGFA on hydrogren radio recombi- 
nation spectral profile linewidtli. The vertical dashed line at 22 
km s~^ represents the thermal linewidth expected for a 10* K gas. 
The horizontal dotted line shows the expected AGFA for random 
orientations. The thick horizontal lines show the weighted mean 
AGFA in the subthermal and suprathermal zones. The grey re- 
gions represent the Icr uncertainties on the means in each zone. 
As linewidths increase beyond 22 km s~^, there is generally better 
alignment between the bubble and magnetic position angles. 

on each average is the weighted dispersion in each re- 
gion. Furthermore, a KS test shows only a 30% chance 
that the AGFA of the subthermal subsample and the 
suprathermal subsample are drawn from the same parent 
population. Therefore, Galactic bubbles with recombina- 
tion linewidths greater than 22 km s~^ are better aligned 
with the average Galactic magnetic field orientation than 
Galactic bubbles with smaller linewidths. 

While the suprathermal sample shows a significant 
(^ 5ct) deviation from the expected mean AGFA for a 
random distribution, the subthermal sample also shows 
a marginal deviation (^ 2a) from that same mean. 

6. DISCUSSION 

The observational evidence presented above shows that 
external magnetic fields are important during the earliest 
phases of the evolution of Hii regions. Having shown 
this, a key following question is what external magnetic 
field strengths are required to generate such effects on the 
morphologies of the bubbles? 

6.1. Magnetic-to-Ram Pressure Ratio 

A first estimate of the requsite magnetic field strengths 
can be found by equating the magnetic and ram pressures 
acting on an expanding ionization front: 



1 



(2) 



where rj is an efficiency factor, pshcii is the density of the 
shell, Mshcii is the expansion speed of the shell, and B is 
the magnetic field strength at the outer boundary of the 
shell. An efficiency factor (77) is included because the ram 
and magnetic energy densities need not be equal for the 
magnetic field to perturb the evolution of the expanding 
front. 

By Bernoulli's principle, a pressure differential in a 
fluid element will give rise to a velocity, u oc ^/dp, so 
that magnetic pressure acts to a decrease the bubble's 
expansion velocity. In the limit of equal ram and mag- 
netic pressures, expansion halts. 

This near equivalence of differential pressure and veloc- 
ity allows the relative pressures of the .Churchwell et ahl 
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(|2006D bubbles to be estimated. The ratio of the major 
and minor axes of each bubble will be equal to the ratio 
of its expansion velocities, 

B~ ~ V~' ^ ' 

Applying Bernoulli's principle and including projection 
effects, the observed bubble axis ratio becomes: 

RvaSi] / Pram '^^^(^)Pmag ^^-^ 

-^niin y Pram Pmag 

where i is the inclination angle of the major axis to 
the plane of the sky (e.g., i = 90° would be a line-of- 
sight along the poles of the bubble, which would ap- 
pear as a circle in projection), Piam is the expansion 
ram pressure, and Pmag is the magnetic pressure affect- 
ing expansion along the minor axes. The predicted ef- 
fects of the magnetic-to-ram pressure ratio and inclina- 
tion angle on the observed bubble axis ratio are shown 
in Figure [71 As the magnetic-to-ram pressure ratio in- 
creases, the projected bubbles are distorted from circles 
(^maj/-Rmin = 1) into ellipscs. Thc inclination angle also 
affects the observed axis ratio. In the limit oi i = 90°, 
all bubbles appear as circles. 
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Fig. 7. — The predicted major-to-minor axis ratio for bubbles 
experiencing different magnetic-to-ram pressure ratios at different 
viewing angles. For i = 90°, the effect is unobservable since the 
anisotropic force acts along the line-of-sight. 

Since inclination will always cause the observed bub- 
ble axis ratio to be smaller than the true axis ratio, bub- 
bles with the largest observed axis ratios (least likely to 
suffer from strong projection effects) can be used to esti- 
mate the typical magnetic-to-ram pressure ratio. Table [5] 
lists the major-to-minor axis ratios and the correspond- 
ing lower limits on the magnetic-to-ram pressure for the 
twelve bubbles with hydrogen recombination linewidths 
greater than 22 km s~^ (those showing the elongation- 
magnetic field correlation). As an example, the largest 
Galactic bubble axis ratio in this sample is 1.38 (bub- 
ble N80). Assuming that this bubble has an inclination 
angle of 0°, then Pmag = 0.47pram (and the coefficient is 
equal to the "q efficiency term). If the inclination is non- 
zero, the magnetic-to-ram pressure ratio is even larger. 
Using the observed bubble axis ratios listed in Table [21 
a lower limit on 77 can be estimated as the unweighted 
mean of these estimates, < 77 >= 0.29. This remains a 
lower limit, because projection effects will cause the ob- 
served axis ratios (and therefore r] values) to be biased 
towards unity, as shown in Figure [71 and so underesti- 
mated ?7true- Uncertainties on bubble major and minor 



TABLE 2 
Magnetic-to-ram 
Pressure Estimates 



Bubble 
Name 
(1) 


Axis 
Ratio 
(2) 


Pressure 
Ratio r) 
(3) 


N23 


1.08 


0.14 


N27 


1.20 


0.31 


N31 


1.30 


0.41 


N50 


1.13 


0.21 


N53 


1.01 


0.02 


N31 


1.22 


0.33 


N66 


1.28 


0.39 


N73 


1.06 


0.12 


N80 


1.38 


0.47 


N96 


1.23 


0.34 


NllO 


1.09 


0.15 


N115 


1.27 


0.38 



axis radii were not reported bv lChurchwell et al.l ()2006D . 
so uncertainties in the axis and pressure ratios were not 
propagated. 

6.2. Magnetic Field Strength Estimates 

A lower limit to the mass of an expanding shell is the 
equivalent mass of the ambient interstellar medium swept 
up by the shell. The shell's density can be estimated as 

PiSM 1 - [(1- < T >)]3 ^ ^ 

where pisM is the average mass density of the interstellar 
medium, and < T > is th e average fractional thic kness 
of the shell, calculated bv lChurchwell et al.l (|2006D . For 
the twelve bubbles with large hydrogen recombination 
linewidths (Av > 22 km s~^), the average shell thick- 
nesses is < T >= 0.24, and Pshcii = l-78pisM- 

Assuming the shell is expanding at ~ 10 km s~^, the 
ion sound speed of a lO'' K gas, into a medium with a den- 
sity equal to that of the diffuse interstellar medium (1.4 
amu cm~'^: [Krumholz fc Matzneril2009[ ) . the surrounding 
magnetic field must be 72 /xG to exert a magnetic pres- 
sure equal to the shell's ram pressure. Using the pressure 
efficiency factor found earlier (77 = 0.29), this decreases 
to 21 ^G. 

This neglects the mass of any denser natal cloud that 
might also be swept up into the expanding shell, which 
would increase the shell mass density. Assume an 04 
main sequence star, with mass of ~ 60 Mq, was at the 
center of a 1 pc radius expa nding H 11 regi on. If the star 
formation efficiency is 30% (|Murravll201lD . there would 
be an additional 200 Mq of natal material in the shell. 
In this case, the magnetic field must be amplified to 2.8 
mG to equal the ram pressure (810 /iG including 77). For 
a 25 Mq 08 main sequence star, this becomes 1.8 mG 
(520 iiG). 

This notion may be supported by OH maser ob- 
servations towards Hii regions. OH masers often 
arise in the shocked neutral gas surrounding ultra- 
compact Hii regions and the magnetic field strength 
can be measured via Zeeman sp litting (jFish fc ReidI 
I2006f) . iDesmurs fc Baudrvl ()1998[ ) observed magnetic 
field strengths of 4 to 6 mG near the ultra-compact H 11 
regions ONI and W51. fish et al., (,2005. ) surveyed 18 



Magnetic Alignment of Galactic Bubbles 



9 



Galactic massive star-forming regions at 1665 MHz with 
the Very Large Baseline Array and observed mG strength 
fields in 184 unique OH masers. Methanol masers have 
also been used to probed mag netic field st r engt hs in high- 
mass star-forming regions. iVlemmingsl (2008) find an 
average magnetic field strength o f 23 mG from 24 b right 
methanol masers. More recentlv. iGreen et al.l ()2012D has 
found mG magnetic fields towards OH masers in high- 
mass star forming regions towards the Carina-Sagittarius 
spiral arm tangent. 

These maser measurements are toward very young H ii 
regions, much smaller and younger than the H ii regions 
interior to Galactic bubbles. However, future high reso- 
lution studies of their morphologies in comparison to the 
local magnetic field orientation (perhaps through high 
resolution polarized dust emission at submillimeter wave- 
lengths) would provide additional insight into the mag- 
netic field properties of H ii region driven Galactic bub- 
bles. 

Detailed calculation of the ram pressure for each bub- 
ble is not possible here. For isolated high mass stars, 
that would require either spectroscopic classification of 
the central ionizing star or direct measurement of the 
mass (i.e., gas column for atomic or molecular species at 
radio wavelengths) and velocity of the expanding shell. 
Spectral classifications are not available for these stars, 
and the assumption of a single, ionizing star may not 
apply if several early-type stars contribute ionizing ra- 
diation. High-resolution spectroscopic observation s with 
airbo rne instruments, such as GREAT (Hevminck et al.l 
I2OII) on SOFIA, may allow more reliable measurements 
of, for example, the 01 (63 /im) traced mass and velocity 
structure around H 11 regions. 

The average magnetic field strength in the solar neigh- 
borhood is around 6 /iG (Beck 2009). In the limit of 
fiux-freezing, and assuming that a magnetic field of this 
strength threads all of the material swept up by a bubble 
shell, the 04 case would need to sweep up a volume hav- 
ing a cross-sectional area equal to 2800 /6 w 470 times its 
current area for the magnetic pressure to equal the ram 
pressure. For a shell of cross-sectional area A = 2-Krdr, 
and assuming that dr = 0.2r, this is an increase in the 
radius of the shell by a factor of 21. In expanding from a 
radius of 0.05 pc to 1 pc, the shell will accumulate enough 
magnetic fiux for the magnetic pressure to equal the ram 
pressure. As discussed earlier, the magnetic pressure is 
typically only a fraction of the ram pressure, decreasing 
the expansion factor necessary to affect the bubble's evo- 
lution. However, relaxing the flux-freezing requirement 
would allow magnetic fiux to diffuse into the bubble and 
(if expanding at sub-alfvenic speeds) possibly away from 
the expanding shell. 

As shown in Figure |6] and discussed in Sec. 15.11 this 
magnetic alignment mechanism must end by the time the 
recombination linewidth falls to 22 km s~^. This roughly 
corresponds to when the inner major axis has expanded 
to approximately 2 pc (see Figure [5]). Magnetic warping 
of bubble shells could also be less important for high mass 
star formation in a dense, extended medium where mass 
loading of the expanding shell causes it to slow before 
appreciable magnetic amplification can occur. 

7. CONCLUSIONS 



The interaction of H 11 region driven Galactic bubbles 
with the Galactic magnetic field has been examined. 
Starlight polarimetry from GPIPS was used to probe 
the large-scale properties of the Galactic disk magnetic 
field and revealed that the projected orientation of the 
Galactic magnetic field is approximately constant ((P- 
GPA) = 85° ± 10°) over the Galactic longitude range 
18° < ^ < 56°, and consistent with previous work show- 
ing the toroidally-dominated nature of the Galactic mag- 
netic field in the disk. Existing hydrogen radio recombi- 
nation line data from HRDS were used to measure phys- 
ical properties of H 11 region driven bubbles and to esti- 
mate their kinematic distances. 

Hii region driven bubbles with suprathermal (> 22 
km s^^) hydrogen recombination linewidths are preferen- 
tially aligned with the average orientation of the Galactic 
magnetic field in the disk. H 11 region driven bubbles with 
subthermal (< 22 km s~^) linewidths are consistent with 
random alignments. Recombination linewidths were also 
shown to be anticorrelated with the physical sizes of bub- 
bles. Since Hii region turbulence should decrease with 
time, recombination linewidths may be used as crude age 
indicators for Hii regions. The spectral types of the stars 
at the centers of the Hii regions were estimated from 
the HRDS kinematic distances and RRL flux densities. 
Spectral type has a weak effect on the observed linewidth 
(but is unable to account for the break at 22 km s^^), 
and showed no correlation with the relative alignment be- 
tween bubble long axes and the mean Galactic magnetic 
field direction (AGFA). 

Th e major-to-minor a xis ra tios of the Galactic bubbles 
from iChurchwell et al.l ()2006[ ) were used to estimate the 
magnetic-to-ram pressure ratio for each bubble. A lower 
limit to this ratio was shown to be 0.29 for the sample of 
magnetically-aligned bubbles (linewidths > 22 km s~^). 

These findings have led to development of a scenario for 
the evolution of the relative alignment between H 11 re- 
gion driven bubbles and the Galactic magnetic field. New 
H II regions are small and characterized by large turbu- 
lent energy densities, showing suprathermal linewidths. 
As H II regions grow and age, their ionization fronts ex- 
pand into magnetized interstellar media, where the shells 
interact with the large-scale Galactic magnetic field. The 
magnetic field is entrained in the shells and becomes am- 
plified. The increasing magnetic field strength inhibits 
expansion perpendicular to the large-scale field orienta- 
tion, causing the bubble to expand preferentially along 
the external magnetic field direction. This creates ellipti- 
cal bubbles preferentially aligned with the magnetic field 
direction. 

As the Hii regions further age, the shells slow be- 
cause of mass loading and the magnetic field weakens 
via magnetic diffusion. Around this time, turbulent en- 
ergy in the hydrogen recombination line emitting region 
dissipates, causing the linewidths to decrease to thermal 
values. Magnetic fields no longer impress a preferential 
orientation on the H 11 region's continued expansion and 
other forces (e.g., local gas density variations) come to 
dominate the morphological evolution of the H 11 region. 
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TABLE 3 

LinUlL.i: I'AHAMETEHS 



Namel Ecc.^ iJ^aj^ B-GPA (P-GPA) 
(arcmin) (deg) (deg) 
(1) (2) (3) (4) (5) 



Outside GPIPS Region 



Nl 


0.47 


0.93 


54 ± 15 


N2 


0.56 


6.96 


108 ± 17 


N3 


0.87 


0.99 


32 ± 3 


N4 


0.53 


1.95 


64 ± 14 


N5 


0.79 


3.71 


39 ± 8 


N6 


0.64 


5.94 


17 ±10 


N7 


0.60 


0.47 


135 ± 3 


N8 


0.00 


0.17 


58 ± 180 


N9 


0.71 


0.59 


133 ± 5 


NIO 


0.80 


1.41 


137 ± 4 


Nil 


0.70 


1.35 


60 ± 4 


N12 


0.50 


4.48 


76 ± 14 


N13 


0.81 


0.57 


75 ± 4 


N14 


0.77 


1.30 


158 ± 5 


N15 


0.75 


1.57 


161 ± 4 


N16 


0.59 


2.42 


35 ± 14 
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TABLE 3 — Continued 



Name^ 


Ecc.i 


p .1 


B-GPA 


(P-GPA) 






(arcmin) 


(deg) 


(deg) 


(1) 


(2) 


(3) 


(4) 


(5) 


N17 


0.68 


0.42 


152 d= 1 




N18 


0.35 


6.61 


64 ± 44 




N19 


0.41 


3.91 


4 ±21 




inside GPIPS Region 




0.66 


1.13 


oz ± O 


CiA A _i_ n 1 
94.4 ih U.i 




0.72 


2.30 


o7 ± O 


oy.o lb U.i 


N22 


0.67 


1.71 


99 ±8 


87.7 ±0.1 


N23 


0.39 


0.39 


19 ±15 


91.4 ±0.1 


N24 


0.47 


7.00 


loo ± 26 


no o 1 o 1 

92. o ± U.I 


Moc: 
IN zo 


0.76 


0.83 


8.8 lb 3 


o v n _L n 1 
O ( .U ± U.I 


N26 


0.63 


0.50 


136 di 3 


89.1 ± 0.1 


N27 


0.56 


1.06 


115 dz 9 


A 1 1 O "1 

94. o IT U.I 


IN zo 


U.40 




VQ 1 1 T 
i6± i 1 


oo. ( zt U. 1 


IN zy 


n 7/1 


9 C9 


97 -U 
Z ( it O 


1 ni _!_ n Q 
lUl.D it U.o 


N30^ 


0.57 


0.99 


55 dz 19 


96.0 ± 0.2 


N31 


0.63 


0.71 


125 d= 6 


A 1 1 n o 

94.1 ± 0.2 


N32 


0.42 


0.37 


180 lb 5 


91.1 ± 0.1 


N33 


0.45 


0.39 


37 ±10 


86.4 ±0.1 


iNo4 


0.60 


1.U7 


37 ± 8 


77. o di U.i 


■\TQr; 
INoO 


0.81 


o.d4 


AO _1_ C 
4z ± O 


yU.4 ih U.i 


N36 


0.75 




OQ _i_ a 
ZO zb D 


oo.y lb U.z 


IN o I 


U.o4 


Z.U 1 


1 -1- /I 
10 it 4 


QQ Q _!_ n 1 
oo.o it U. 1 


IN OO 


fl ^9 
U. OZ 


U. Oo 


ID ( III 10 


on A 4- fl 1 

oU.4 31 U. 1 


IN oy 


U.DO 


1 .yo 


/ o III y 


89 fl -1- fl 9 
OZ.U III U.Z 




u. / u 


1 9*^ 
1 .ZO 


1 1 m o 


8Q A -1- fl 1 
oo.'y: III U. 1 




fl 9^ 

u.zo 


fl 49 


12-1- 7fl 
lo III 1 U 


80 A -1- fl 1 

OU.'rt IC U.I 


N42 


0.76 


0.67 


150 ±6 


83.9 ± 0.1 


N43 


0.58 


0.61 


1.3 ±9 


89.1 ± 0.1 


KIAA 


U. 00 


1 1 fl 
i . lU 


1 fic; _|_ Q 
IDO It y 


Q1 A -1- fl 1 

y i.4t It U.I 


N45*^ 


0.78 


1.53 


50 ±15 


riri o 1 o o 
))().() ± U.Z 


IN 40 


U.oo 


1 T O 

i.iy 


O Q -1- 91 

y.o it zi 


Dy.D it U.Z 




U. 1 O 


9 Af\ 

Z.'iU 


04 III D 


luu.z m U.Z 


N48 


u.oo 


u.o / 


142 ± 18 


1 1 n 9 _1_ fl 9 

XJ.U.Z HI U.Z 


iN'iy 


U.Do 


1 ^9 

l.OZ 


on -1- Q 
yu m y 


1 0A Q -1- O ^ 
lu^.y m u.o 


IN ou 


U.4D 


i .oy 


9zt -1-17 

z^y: mi/ 


1 1 f; A -1- O *? 

iio.^y: It U.O 


IN Oi 


U.ZO 


i. ( O 


4D it D4 


01 Q J- n A 

yi.o it U.4 


N52 


0.79 


2.22 


48 ±5 


76.5 ± 0.2 


N53 


0.15 


0.74 


134 ± 133 


83.7 ±0.2 


IN o4 


0.70 


1.88 


80 ± 16 


ft v /I _L n o 
D ( .4 ± U.Z 


N55 


0.67 


0.77 


113 ± 6 


nn o _L n a 

yy.z ± U.4 


N56 


0.19 


0.99 


T ,1 1 1 no 

74 ± lUz 


on 1 loo 

yu.i ± U.Z 


iNo / 


0.57 


0.31 


75 ± 7 


O /I Q _L n o 
o4.o lb U.Z 


N58 


0.62 


0.17 


178 ± 1 


on r' loo 

yu.D ± U.Z 


N59 


0.68 


7.03 


23 ± 9 


87.1 ±0.1 


iNoU 


U.DO 


U.D4 


Loo ± 4 


7n Q _i_ n Q 
(U.o lb U.o 


iMOl 


u.ou 


"i 98 


1 98 -1- fi9 
IZo ZL DZ 


7Q 9 -1- O ^ 
(y.z m u.o 


N62 


0.38 


1.41 


36 ±26 


75.6 ± 0.2 


N632 


0.68 


11.89 


25± 16 


90.5 ±0.1 


N64 


0.69 


5.08 


42 ± 9 


OO 1 loo 

93.1 ± 0.2 


N65 


0.49 


2.02 


10 ± 16 


81.2 ± 0.2 


N66 


0.62 


0.46 


90 ± 7 


94.3 ± 0.2 


N67 


0.38 


0.75 


95 ± 13 


103.0 ± 0.1 


N68 


0.72 


5.17 


9.0 ± 10 


94.2 ±0.1 


N69 


0.00 


8.42 


6.3 ± 180 


93.5 ±0.1 


N70 


0.46 


0.50 


38 ± 17 




N71 


0.53 


7.08 


74± 16 


95.5 ± 0.1 


N72 


0.36 


0.92 


41 ± 26 


92.0 ± 0.2 


N73 


0.34 


0.77 


41 ± 24 


84.1 ±0.2 


N74 


0.59 


1.27 


52 ±7 


90.7 ±0.3 


N75 


0.57 


0.62 


172 ±6 


88.9 ±0.3 


N76 


0.65 


3.73 


77 ±10 


96.1 ±0.2 


N77 


0.33 


1.16 


74 ±31 


68.7 ±0.1 


N782 


0.68 


0.31 


135 ± 16 


71.1 ±0.1 


N79 


0.64 


1.27 


169 ±9 


65.4 ±0.2 


N80 


0.69 


1.90 


62 ±9 


72.1 ± 0.2 


N81 


0.70 


9.24 


73 ±7 


64.8 ±0.1 


N82 


0.13 


1.47 


119 ± 180 


77.3 ± 0.3 
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Eec.i 




B-GPA 


(P-GPA) 






(arciilin) 






(1) 


(2) 


(3) 


(4) 


(5) 


NSS^ 


0.65 


0.31 


155 ± 17 


73.1 ± 0.2 


N84 


0.69 


1.19 


3.1 ±6 


77.0 ±0.3 


N85 


0.57 


0.46 


104 ±3 


85.3 ± 0.3 


N86 


0.46 


0.26 


50 ± 2 


85.2 ± 0.3 


N87 


0.61 


0.31 


32 ±6 


87.6 ± 0.3 


N88 


0.30 


1.37 


176 ± 32 


94.7 ±0.4 


N89 


0.47 


0.92 


154 ± 15 




N90 


0.13 


1.53 


122 ± 180 




N91 


0.74 


5.19 


40 ± 8 


107.1 ± 0.3 


N92 


0.73 


2.03 


17 ±7 




N93 


0.68 


0.73 


9 ± 8 




N94 


0.74 


3.85 


96 ±6 


85.3 ± 0.2 


N95 


0.60 


2.13 


125 ±9 


87.2 ± 0.5 


N96 


0.59 


0.40 


11 ±3 


97.0 ±0.2 


N97 


0.67 


4.23 


95 ±8 


86.0 ±0.6 


N98 


0.11 


1.37 


43 ± 180 


97.2 ± 0.3 


N99 


0.64 


4.20 


177 ± 11 


95.4 ±0.3 


NlOO 


0.63 


4.85 


54± 11 


82.1 ± 0.2 


NlOl^ 


0.52 


0.99 


125 ± 21 


81.9 ± 0.2 


N102 


0.12 


1.82 


171 ± 180 


82.2 ± 0.3 


N103 


0.63 


0.63 


91 ±9 


82.0 ±0.3 


N104 


0.64 


0.77 


34 ±8 


81.0 ±0.2 


N105 


0.75 


0.87 


103 ±5 


79.8 ± 0.2 


N106 


0.65 


0.39 


36 ± 10 


84.8 ± 0.2 


N107 


0.62 


11.59 


110 ± 12 


81.1 ± 0.1 


N108 


0.84 


4.39 


137 ±4 


77.0 ±0.2 


N109 


0.61 


14.57 


109 ± 12 


78.7 ±0.1 


NllO 


0.38 


0.55 


122 ± 19 


71.6 ± 0.1 


Nlll 


0.59 


1.00 


44 ± 8 


71.4 ±0.1 


N112 


0.74 


0.28 


143 ±7 


71.6 ± 0.1 


N113 


0.80 


0.57 


38 ±3 


72.4 ±0.1 


N114 


0.20 


1.35 


72 ± 75 


71.9 ±0.1 


N115 


0.62 


3.20 


80± 12 


81.3 ± 0.2 


N116 


0.61 


0.79 


118 ±5 


77.8 ± 0.2 


N117 


0.14 


1.41 


48 ± 180 


77.4 ±0.2 


N118 


0.58 


0.32 


122 ±6 


74.1 ± 0.2 


N119 


0.48 


4.31 


178 ± 20 


77.6 ± 0.2 


N120 


0.33 


1.26 


99 ± 32 


72.7 ±0.3 


N121 


0.55 


0.49 


144 ±9 


70.2 ± 0.4 


N122 


0.72 


0.53 


177 ±3 


79.1 ± 0.4 


Outside GPIPS Region 


N123 


0.73 


1.47 


106 ±6 




N124 


0.35 


1.60 


161 ±33 




N125 


0.44 


0.74 


142 ± 19 




N126 


0.61 


1.82 


69 ±9 




N127 


0.20 


2.94 


24 ± 178 




N128 


0.70 


3.38 


60 ±8 




N129 


0.70 


3.27 


13 ± 17 




N130 


0.47 


0.98 


169 ± 13 




N131 


0.47 


6.18 


77 ±21 




N132 


0.36 


0.21 


86 ± 7 




N133 


0.78 


1.83 


73 ± 5 




N134 


0.83 


0.58 


79 ± 3 





1 Data taken from lChurchwell et al] HWOB) . 



^ Bubble fitting failed, CPA measured as described in text. 



